arXiv:hep-ph/9706477vl 24 Jun 1997 


PRINCE TON/HEP/97-5 
|hep-ph/9706477 

An Overview of D°D° Mixing Search Techniques * 


Tiehui (Ted) Liu 

Department of Physics, Princeton University, Princeton, NJ 08544 


Abstract 

The search for D°D° mixing may carry a large discovery potential for new 
physics since the D°D° mixing rate is expected to be small in the Standard 
Model. This paper gives a brief review of the experimental techniques which can 
be used to search for D°D° mixing. 


1 Introduction 

Much of the enthusiasm for searching for D°D° mixing stems from the belief that the 
search carries a large discovery potential for New Physics, since the mixing rate Rmixing = 
B(D° —* D° —> /)/ B(D° —> /) is expected to be small in the Standard Model (range from 
10 -7 to 10 -3 , see references in elsewhere [p] ||. ||). One can characterize D°D° mixing in 
terms of two dimensionless variables: x = 8m/y+ and y = y-/y+, where the quantities y± 
and 8m are defined by y± = (71 ± 72)/2 and 5m = m 2 — m\ with m,;, 7 i (i = 1 , 2 ) being 
the masses and decay rates of the two CP (even and odd) eigenstates. Assuming a small 
mixing, namely, 8m, y- <C 7 + or x, y <C 1, we have R m bdng = ( x 2 + y 2 )/ 2 . An overview 
of the current experimental status and future prospects can be found elsewhere |jj. This 
paper only gives a brief review of the experimental techniques which can be used to search 
for D°D° mixing. 


2 The Techniques 

The techniques which can be used to search for mixing can be roughly divided into two 
classes: hadronic and semi-leptonic. Each method has advantages and limitations, which 
are described below. 
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2.1 Hadronic Method 

The hadronic method is to search for the D° decays D° —> I\ + n~ (X) Jjj. These decays 
can occur either through D 0 D° mixing followed by Cabibbo favored decay D° —> D° —> 
A' + 7 T _ (X), or through DCSD (Doubly Cabibbo Suppressed Decay) D° —> A' + 7 r“(X). This 
means that the major complication for this method is the need to distinguish between DCSD 
and mixing. The hadronic method can therefore be classified according to how DCSD and 
mixing are distinguished. In principle, there are at least three different ways to distinguish 
between DCSD and mixing candidates experimentally: (A) use the difference in the decay 
time-dependence; (B) use the possible difference in the resonant substructure between DCSD 
and mixing events in D° —> , etc. modes; (C) use the quantum 

statistics of the production and the decay processes; (D) there is also a possibility of using the 
CP eigen states, such as D° —> K + K~, Tr + n~ (Singly Cabibbo Suppressed Decays - SCSD), 
to study mixing caused by the decay rate difference: y = 7-/7+ = (72 — 7 i)/( 7 i +72)- Note 
that DCSD decays are interesting in their own right, since they can be used to extract CKM 
angle 7 in B~ —> D°(D°)K~ modes ||. 


2.1.1 Method A -use the difference in the time-dependence of the decay 

This method is to measure the decay time of the D° —> K + tt~ decay. Here the D° 
tagging is usually done by using the decay chain D* + —> D°ir+ followed by D° —> K + tt~. 
A pure D° state generated at t = 0 decays to the A' + 7 t“ state either by D°D° mixing or by 
DCSD, and the two amplitudes may interfere. 

Following the notation in Ref.l, and define a(f) = Amp(D° —■> /), a(f) = Amp(D° —> 
/) with p(f) = a(f)/a(f); a{f) = Amp(D° -> /), a (/) = Amp(D° -> /) with p(f) = 
a(f)/a(f), and 77 = |p(/) (fj = pP(f) ), for 6mt,6'yt « 1 and small \p\ ( in the case of 
/ = AT+7T - ), we have [0 

I( \D° phys (t) >- / ) = Hf) | 2 |^| 2 | \{i8m+ ^ 7 ) t + y | 2 e~^. (1) 


\D° phys (t) >- / ) = \a(f)\ 2 |^| 2 | ±(i5m+ l -5y) t + fj | 2 e~^. (2) 


Note the difference between Equation |l] and is the indication of CP violation. Next let 
us write the equations in a somewhat different form, which is more convenient for discussion 
here. Note that j |( iSm + ^ 7 ) t + 77 | 2 can be written in the form: | ^(ix + y) t + 1 7 | 2 
where now the time t is measured in unit of average D° lifetime (I/ 7 +). Recall that we 
have Rmixing = ( x 2 + y 2 )/ 2 and define Rdcsd = \f)\ 2 , the term above can be written in the 
form | \/Rmixing/2 t + VR-dcsd ^ | 2 where <f> = Arg(ix + y) —Arg(rf) is an unknown phase 
between the two amplitudes. 


Equation jl] and now become: 

I(\D° phys (t) >/) = |d(/)| 2 |^| 2 x 


Rdcsd T y - RmixingR DCSD t COS 0 7 ^ Rmixing 


( 3 ) 


and 

J(|D° phys (f) >—»• /) = |a(/)| 2 |^| 2 x 
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(4) 


Rdcsd + v 2R, 


s Rdcsd t cos (j) ■ 


~ Rmixing t 


where Rdcsd = \v\ 2 and 4> = Arg{ix + y) — Arg(fj). Note (f> is different from cj> because 
Arg(fj) = Arg(q/p) + Arg(p(f)) where Arg(q/p) = —Arg(p/q ) (and one can define q/p = 
e - 2 i<j> M j n the usual way). For convenience, let us ignore the CP violation in the decay 
amplitude. Thus the difference between the interference terms (that is, the interference 
phase <j> and (j)) in Equation and ^ would be the indication of CP violation. 

Assuming CP conservation, Equation ^ and [| simply become Jg| : 


IG-DphysW —»■ K + tt )(t) oc (Rdcsd + v^RmixingR dcsd t cos<j) + 2 -^mixing^ *■ (5) 

The signature of mixing is therefore a deviation from a perfect exponential time dis¬ 
tribution with the slope of 7 + . A small mixing signature could show up in the interference 
term at lower decay times. The importance of the interference term has been discussed in 
detail elsewhere |], |l|. The importance of possible large CP violation effect (due to New 
Physics) has been discussed by Wolfenstein [jTj and others Q. It has also recently been ar¬ 
gued by Browder and Pakvasa || that it may be possible to calculate the phase Arg(p(f )) 
due to final state interaction. 

It is interesting to take a look at the time integrated effect. Let us define 

p _ J 0 °°n\n° hys (t)>^f)dt 
f 0 °°I(\DO phys (t) -/),/,• 

and we have |ij 

R — Rdcsd T \/2H tm jxingR dcsd cos (j) T Rmixing* (7) 

In the special case when | cos(^| = 1, Equation can be written in the form 

( Rdcsd — R ) 2 + ( Rmixing — R ) 2 = R 2 . (8) 


2.1.2 Method B -use the difference in the resonant substructure between DCSD and mixing 
events in multi-body decays 

The idea |(>|| is to use the wrong sign decay D* + —> D° 7r+ followed by D° —> A ,+ 7T _ 7r°, 
I\ + n~n + 7T~ , etc., and use the possible differences of the resonant substructure between 
mixing and DCSD to study mixing. In principle, one can use the difference between the 
resonant substructure for DCSD and mixing events to distinguish mixing from DCSD. For 
instance, combined with method A, one can perform a multi-dimensional fit to the data by 
using the information on AM, M(D°) 1 proper decay time t and the yield density on Dalitz 
plot n w (p,t). The extra information on the resonant substructure will, in principle, put a 
better constraint on mixing. Of course, one needs large amount of clean data to do this 
in the future. Because of this, for current experiments this method is more likely to be a 
complication rather than a better method. In general, we cannot treat multi-body decays 
exactly the same way as D° —> AT + 7r - when the resonant substructure is ignored. One can 
find detail discussions about this elsewhere [|Ijj. 
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2.1.3 Method C —use quantum statistics of the production and decay processes 


This method is to search for dual identical two-body hadronic decays in e + e - —> T" —> 
D°D°, such as (K~tt + )(K~tt + ), as was first suggested by Yamamoto in his Ph.D thesis [|Tot . 
The idea is that when D°D° pairs are generated in a state of odd orbital angular momentum 
(such as T"), the DCSD contribution to identical two-body pseudo-scalar-vector (D —> PV) 
and pseudo-scalar-pseudo-scalar (D —> PP) hadronic decays (such as (K~n + )(K 
cancels out, leaving only the contribution of mixing @0 in the wrong sign versus right 
sign ratio: 


P _N (K tt+,A 7t+)+N(A+7t ,K+n ) ( 2 o . /0 p 

K “ N(A-tt+, K+n~) + N(A+tt-, K-n+) ~ [X + V ~ Kmixing 


(9) 


One common misunderstanding about this method is that DCSD is forbidden in this case. 
This is not true, DCSD is allowed with mixing and contributes to both wrong sign and right 
sign (but the contribution cancels out in the ratio). One can find the essence of Yamamoto’s 
original calculation for the (A' _ 7 t + )(A' _ 7 t + ) case elsewhere 0]. Detailed calculations for this 
kind of method can be found in a recent paper |l 2 |. 


2.1.4 Method D - use CP eigen final states (such as D° —> K + K , 7 r + 7 r ) to measure the 
decay rate difference y 

This is because (assuming CP conservation) those decays occur only through the CP 
even eigenstate, which means the decay time distribution is a perfect exponential with the 
slope of 7 i. Therefore, one can use those modes to measure 71 . The mixing signature is not 
a deviation from a perfect exponential (assuming CP conservation), but rather a deviation 
of the slope from (71 + Since 7 + = (71 + 72)/2 can be measured by using the 

D° —> K~n + decay time distribution, one can then derive y = 7 -/ 7 + = (72 — 7 i)/( 7 i + 72 )- 
Observation of a non-zero y would demonstrate mixing caused by the decay rate difference 
(Rmixing = (x 2 + y 2 )/ 2). Note that there is no need to tag the D°, since we only need to 
determine the slope. 


2.2 Semi-Leptonic Method 

The semi-lcptonic method is to search for D° —> D° —> Xl~n decays, where there is 
no DCSD involved. However, it usually (not always) suffers from a large background due 
to the missing neutrino. In addition, the need to understand the large background often 
introduces model dependence. New ideas are needed in order to improve the sensitivity 
significantly. One can find more discussions elsewhere [jjj |l3| . 


3 Comments On Experimental Issues 

Generally speaking, vertexing and tagging are the two critical elements to the study of 
mixing and CP violation. For BdBd and B S B S: tagging has always been one of the major 
issues and has been proven to be difficult. In D°D° case, tagging seems to be easier as one 
has the clear advantage of using the decay chain D* + —► D°ir+ (see below), however, the 
much smaller D°D° mixing rate requires much cleaner tagging. 

As mentioned before, the D° tagging is usually done by using the decay chain D* + —> 
D°tt+ followed by D° —> A' + 7 r“. The 7 r+ from D* + has a soft momentum spectrum and 
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is referred to as “the slow pion”. The charge of the slow pion is correlated with the charm 
quantum number of the D° meson and thus can be used to tag whether a D° or D° meson 
was produced. The idea is to search for wrong sign D* + decays, where the slow pion has 
the same charge as the kaon arising from the D° decay. The right sign signal D* + —> D°tt+ 
followed by D° —> K~ir + can be used to provide a model-independent normalization for 
the mixing measurement. The small Q value of the D* + decay results in a very good mass 
resolution in the mass difference AM = M(D* + ) — M(D° ) — M( 7r+) and allows a D* + 
signal to obtained with very low background. In general, mixing search requires at least: 

(1) excellent vertexing capabilities, at least good enough to see the interference structure; 

(2) low background around the primary vertex. 

The vertexing capabilities at e + e _ experiments (C/a ~ 3) for CLEO III and asym¬ 
metric B factories at SLAC and KEK may be sufficient for a mixing search. The extra 
path length due to the Lorentz boost, together with the use of silicon detectors for high 
resolution position measurements, have given the fixed target experiments an advantage in 
vertex resolution (typically C/a ~ 8 — 10) over e + e _ experiments. 

The background level around the primary vertex could be an important issue. Low 
background around primary vertex means that one does not suffer much from random slow 
pion background. Define B as the number of wrong sign background events, S the number 
of signal events and Q the D* — D mass difference, one has [[b|: 


B 

~S 


, 1 dB 


( 10 ) 


where oq is the mass difference resolution. Note the background density, is a 

characteristic of the fragmentation process. At CLEO II, the background density is about 
0.001 per MeV and ctq ~ 0.7 MeV which means that B/S ~ 10 -3 . For CLEO III and BF, 
oq can be reduced down to 0.3 MeV with the silicon vertex detector. The low background 
around the primary vertex at e + e~ experiments is a certain advantage. One disadvantage at 
fixed target experiments is the higher background around the primary vertex (higher B/S). 

While vertexing capabilities will likely be improved in the future, the background den¬ 
sity seems difficult to improve. However, there are some ideas Jg| to improve the background 
density for the future experiments. For example, for asymmetric B factory or Z factory, fixed 
target experiments (such as HEAR-B) and especially hadron machines (TEVATRON, LHC), 
it maybe possible to use B° —> D* + l~v, where the primary ( D* + decay) vertex can be de¬ 
termined by the l~ together with the slow pion coming from the D* + . In this case, the 
background level around the primary vertex is intrinsically low since the D* (or B) decay 
vertex is no longer the event vertex. Here the background density, is no longer a 

characteristic of the fragmentation process. In fact, CDF jbij has recently obtained rather 
clean right sign signal for B° —> D* + l~v followed by the decay chain D* + —■> D° 7r+ and 
D° — > K~tt + (~ 1000 signal events from llOpb" 1 data). 

It maybe also possible to use hadronic B decays (assuming one can trigger on these 
events), such as B° —► D* + n~ (BR ~ 0.3%), B° —> D* + 7r _ 7r° (BR ~ 1.5%), B° —> 
D* + n~Tr + n~ (BR ~ 1.2%) and B° —> D* + 7r _ 7r + 7r _ 7r° (BR ~ 3.4%). In addition, charged 
B decays such as B~ —> (BR ~ 0.2%), B~ —> (BR ~ 1.5%) can 

be also used. There should be enough tracks which can be used to determine the B (or 
D* + ) decay vertex. Here not only is the background level around the D* + primary vertex 
intrinsically low, but also the backgrounds can be further reduced by the B mass constraint 
as one can fully reconstruct the B decays (unlike in the B° — > D* + l~v case). Note that 
this is similar to the idea |L5) of using B c —> B s tt + (p + ) to tag B s . One major difference 
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is that B c production cross section is unknown and could be very low, while one expects 
more than 10 10 B produced at Tevatron and LHC which can be potentially used for D°D° 
mixing search. 


4 Summary 

We have briefly discussed some possible techniques which can be used for mixing 
searches in the future. They can be roughly divided into two classes: hadronic and semi- 
leptonic. Each method has advantages and limitations. In the case of semi-leptonic method, 
the advantage is that it is theoretically clean as there is no DCSD involved while the limita¬ 
tion is that very often it is limited by large background due to the missing neutrino. In the 
case of D° —» K + ir~ (X), the major complication is the presence of DCSD. However, this 
very complication can be turned to advantage since the potentially small mixing signature 
could show up in the interference term (at lower decay times). The design of future exper¬ 
iments should focus on improving the vertexing capabilities and reducing the background 
level around the primary vertex, in order to fully take advantage of having the possible 
DCSD and mixing interference. In addition, we have learned that the very complication 
due to the possible differences between the resonant substructure in many DCSD and mixing 
decay modes D° —> K + n~X could also, in principle, be turned to advantage by providing 
additional information. 

In the case of D° —> K + n~(X ) and D° —> X + l ~, we are only measuring Rmbdng = 
(. x 2 + y 2 )/ 2. Since many extensions of the Standard Model predict large x = 8m/ 7 +, 
and we expect New Physics does not affect the decays in a significant way thus does not 
contribute to y, it is important to measure x and y separately. Fortunately, decays such as 
D° —> K + K ~, 7 r + 7 r“, can provide us information on y. Observation of a non-zero y would 
demonstrate mixing caused by the decay rate difference. This, together with the information 
on Rmixing obtained from other methods, we can in effect measure x. It is worth to point 
out here that x and y are expected to be at the same level within the Standard Model, 
however we do not know for sure exactly at what level since theoretical calculations for 
the long distance contribution are still plagued by large uncertainties. Therefore, it is very 
important to measure y in order to understand the size of x within the Standard Model, 
so that when D°D° mixing is finally observed experimentally, we will know whether we are 
seeing the Standard Model physics or new physics beyond the Standard Model. 


Acknowledgement 

I would like to express special thanks to Lincoln Wolfenstein for many useful and stim¬ 
ulating discussions at this conference and Honolulu CP conference. This work is supported 
by Robert H. Dicke Fellowship Foundation at Princeton University. 


References 

References 

[1] T. Liu, “An Overview of D°D° Mixing Search Techniques: Current Status and Future 
Prospects,” in Proceedings of the Workshop on the Tau/Charm Factory, page 447-479, 


6 



Argonne National Laboratory, Argonne, IL, June 21-23, 1995. Princeton/HEP/95-6, 
|hep-ph/9508415| , and references therein. 

[2] S. Pakvasa, these proceedings. 

[3] E. Golowich, “Theory: Rare Decays, Mixing, CPV in D Mesons,” in Proceedings of the 
B Physics and CP Violation Conference, University of Hawaii, Honolulu, March 24-27, 
1997. 

[4] We discuss D° decays explicitly in the text, its charge conjugate decays are also implied 
throughout the text unless otherwise stated. 

[5] D. Atwood, I. Dunietz, A. Soni, Phys. Rev. Lett. 78, 3257 (1997). 

[6] T. Liu, “The D°D° Mixing Search - Current Status and Future Prospects,” in Proceed¬ 
ings of the Charm 2000 Workshop, FERMILAB-C onf-94/190, page 375-394, Fermilab, 
June 7-9, 1994. Harvard preprint HUTP-94/E021, |hep-ph 940833C . 

[7] L. Wolfenstein, Phys. Rev. Lett. 75, 2460 (1995). 

[8] G. Blaylock, A. Seiden and Y. Nir, Phys. Lett. B 355, 555 (1995). 

[9] T. E. Browder and S. Pakvasa, Phys. Lett. B 383, 475 (1996). 

[10] H. Yamamoto, Ph.D thesis, CALT-68-1318 (1985). 

[11] I. Bigi and A. I. Sanda, Phys. Lett. B 171, 320 (1986). 

[12] Z. Xing, Phys. Rev. D 55, 196 (1997). 

[13] R. J. Morrison, “Charm 2000 Workshop Summary,” in Proceedings of the Charm 2000 
Workshop, FERMILAB-Conf-94/190, page 313-322, Fermilab, June 7-9, 1994. 

[14] The CDF II Collaboration, “The CDF II Detector Technical Design Report,” 
FERMILAB-Pub-96/390-E, October, 1996. 

[15] T. Liu and S. Pakvasa, “Rare Decay Searches at CO,” in Proceedings of the Workshop 
on Heavy Quark Physics at C-Zero, Fermilab, Dec. 4-6, 1996. 


7 







